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Abstract 

A lattice QCD calculation of the B -^ Dlv decay form factors is presented. 
We obtain the value of the form factor h+{w) at the zero-recoil limit w = 1 
with high precision by considering a ratio of correlation functions in which 
the bulk of the uncertainties cancels. The other form factor h-{w) is cal- 
culated, for small recoil momenta, from a similar ratio. In both cases, the 
heavy quark mass dependence is observed through direct calculations with 
several combinations of initial and final heavy quark masses. Our results 
are h+{l) = 1.007(6)(2)(3) and /i_(l) = -0.107(28)(04)(^[J). For both the 
first error is statistical, the second stems from the uncertainty in adjusting 
the heavy quark masses, and the last from omitted radiative corrections. 
Combining these results, we obtain a precise determination of the physical 
combination JFq^/)(1) = 1.058(^y), where the mentioned systematic errors 
are added in quadrature. The dependence on lattice spacing and the effect 
of quenching are not yet included, but with our method they should be a 
fraction of J-'b^d — 1- 



1 Introduction 

The precise determination of the Cabibbo-Kobayashi-Maskawa (CKM) ma- 
trix element Vet is a crucial step for B physics to pursue phenomena beyond 
the Standard Model. In particular, the precision achieved in determining 
the apex of the unitarity triangle may be limited by \Vcb\, even with future 
high-statistics experiments. The current determination of \Vcb\ [0] is made 
through inclusive p|,|^ and exclusive [^|^ B decays. 

The heavy quark expansion offers a method to evaluate the hadronic 
transition amplitude in a systematic way. In particular, at the kinematic end 
point the exclusive B ^ D* matrix element is normalized in the infinite heavy 
quark mass limit, and the correction of order l/mq vanishes as a consequence 
of Luke's theorem |^. It is thus possible to achieve an accuracy on \Vcb\ of a 
few percent. Calculations of the ^/ttlq (and higher order) deviations from the 
heavy quark limit have previously been attempted with the non-relativistic 
quark model and with QCD sum rules. 

Lattice QCD has the potential to calculate exclusive transition matrix 
elements from first principles. The shapes of the B -^ D^*'>lu decay form 
factors have already been calculated successfully with propagating |0-0, 
static ]T(]|^r5[, and non-relativistic |]n[ heavy quarks. On the other hand, 



a precise determination of the absolute normalization of the form factors has 
not been achieved. This paper fills that gap for the decay B —>■ DW. 

Previous lattice calculations were unable to obtain the normalization of 
the form factors for various reasons. First, the statistical precision of the 
three point function (DVf^B'^), which is calculated by Monte Carlo integra- 
tion, has not been enough. Second, perturbative matching between the lat- 
tice and the continuum currents has been a large source of uncertainty. Since 
the local vector current defined on the lattice is not a conserved current at 
finite lattice spacing a, the matching factor is not normalized even in the 
limit of degenerate quarks. Although one-loop perturbation theory works 
significantly better with tadpole improvement ITS , the two-loop contribution 



remains significant (a^ ~ 5 %). Last, the systematic error associated with 
the large heavy quark mass must be understood. Previous work with Wilson 
quarks 0^, for which the discretization error was as large as O^arriQ), could 
not address the l/mq dependence in a systematic way when mq ^ 1/a. 
In this paper we present a lattice QCD calculation of the B -^ DW 



decay form factor. For the heavy quark we use an improved action |16 



for Wilson fermions, reinterpreted in a way mindful of heavy-quark symme- 



try |T^. Discretization errors proportional to powers of aniq do not exist in 
this approach. Instead, discretization errors proportional to powers of oAqcd 
remain, although they are intertwined with the l/rnq expansion. The first 
extensive application of this approach to heavy-light systems was the calcu- 
lation |l^,[l9| of the heavy-light decay constants, such as fs and fo- There 
the lattice spacing dependence was studied from direct calculations at sev- 
eral lattice spacings, and a very small a dependence was observed. The third 
difficulty mentioned above is, thus, no longer a problem. 

To obtain better precision on the semi-leptonic form factors, we introduce 
ratios of three-point correlation functions. The bulk of statistical fluctuations 
from the Monte Carlo integration cancels between numerator and denomina- 
tor. Furthermore, the ratios are, by construction, identically one in both the 
degenerate-mass limit and the heavy-quark-symmetry limit. Consequently, 
statistical and all systematic errors, as well as the signal, are proportional 
to the deviation from one. The first and second difficulties given above are, 
thus, also essentially cured. 

The ratio of correlation functions for the calculation of h^{l) corresponds 
to the ratio of matrix elements, 

{D\c^ob\B){B\hoc\D) _ 

{D\-c^oc\D){B\hob\B) l^+^'^l ' ^'^ 

in which all external states are at rest. The denominator may be considered 
as a normalization condition of the heavy-to-heavy vector current, since the 
vector current q'~ffj,q with degenerate quark masses is conserved in the con- 
tinuum limit, and its matrix element is, therefore, normalized. As a result 
the perturbative matching between the lattice and continuum currents gives 
only a small correction to |/i+(l)|. 

For the calculation of h_{w) we define another ratio, corresponding to 
matrix elements 

{D\ciih\B) {D\c^,c\D) ^ ^ _ h_{w) 

{D\c^^h\B) {D\c^,c\D) h+{wy ^^ 

where equality holds when the final-state D meson has small spatial mo- 
mentum. By construction, the ratio produces a value of /i_ that vanishes 
when the b quark has the same mass as the c quark, as required by current 
conservation. 



This method does not work as it stands for the B -^ D*W decay form 
factors. The axial vector current mediates this decay, and it is neither con- 
served nor normahzed. We will deal separately with this case in another 
paper. 

This paper is organized as follows. Section ^ contains a general discussion 
of form factors for the exclusive decay B —>■ DW. Sections ^ and ^ discuss 
heavy quark effective theory and the l/mg expansion in the continuum and 
with the lattice action used here. Section ^ contains details of the numerical 
calculations. Sections ^-^ present our results. Sections ^ and |^ discuss the 
form factor h^ and its mass dependence. Sections |^ and ^ do likewise for /i_. 
We compare the results from the fits of the mass dependence to corresponding 
results from QCD sum rules in Sec. ITUl The values of /i+(l) and h_{l) at 



the physical quark masses are combined in Sec. into a result for the form 
factor J^b^d{^), which with experimental data determines \Vcb\- We give our 



conclusions in Sec. 12 



2 B ^ Dlu form factors 

The decay amplitude for B ^ Dli? is parametrized with two form factors 
h^{w) and h^{w) as 

(D(p')|V^|5(p)) = ,/TK^^ [hl--''{w){v + v')^ + h''_-'>{w){v - v')^] , (3) 

where v and v' are the velocities of the B and D mesons, respectively, and 
w = V ■ v'. The square of the momentum transferred to the leptons is then 
g^ = m^ + m|) — 2mBm£)W. We denote by the symbol V^ the physical 
vector current, to distinguish it from currents in heavy quark effective theory 
(HQET) and in lattice QCD. 

The differential decay rate reads 

aw ASiT-^ 

with 

:Fb^d{w) = hl-^{w) - ^^~^^ /^^-^H. (5) 

■niB + mD 

At zero recoil {v' = v, so w = 1) one expects Tb^d{X) to be close to one, 
because of heavy quark symmetry. From (HI) a determination of \Vcb\ consists 



of the following three steps: measure iKbHJ^B^Dl'W^)! in an experiment, ex- 
trapolate it to the zero-recoil limit assuming some functional form, and use 
the theoretical input of J-'b-,d{^)- 

In this paper we report on a new calculation of JFb^d(I) with lattice 
QCD, which is model independent, at least in principle.^ The present cal- 
culation includes the leading corrections to the heavy-quark limit: radiative 
corrections to the static limit of h^~^^{l), the l/rnQ contribution of h^^^{l), 
and the l/mg contributions of /i^^^(l). Radiative corrections of order as 
to ^-(1) are not yet available, but these and further corrections, of order 
ttg, as/rriQ, etc., could be included in future applications of our numerical 
technique, once the needed perturbative results become available. 

An obvious disadvantage in using the B -^ DW decay mode is that the 
branching fraction is much smaller than the B -^ D*W mode. Another, but 
not less important, shortcoming is that the phase-space suppression factor 
{w'^ — 1)^/^ makes the extrapolation of the experimental data to w = 1 
more difficult than for B -^ D*W, where the corresponding factor is [w"^ — 
1)^/^. Nevertheless, the experimental result of the CLEO Collaboration pO[] 
shows that the above method certainly works, even with current statistics. 
That means that future improvement of the statistics will allow a much 
better determination of \Vcb\, providing an important cross check against 
other methods. 



3 HQET and the l/mg expansion 

Many important theoretical results have been obtained for the form factors 
with HQET. The Lagrangian of HQET uses fields of infinitely heavy quarks, 
so that heavy quark symmetries are manifest. The effects of finite quark mass 
are included through the l/mq expansion and through radiative corrections. 
For example, at zero recoil the form factor h+ is given by 



h+{l) =riv 






(6) 



where rjv represents a matching factor relating the vector current in (^) to 
the current in HQET [^. The absence of the 0(1/?tiq) term in (^ is a 

"'^Our calculations are done in the quenched approximation, for example, but this is a 
removable uncertainty and not a permanent limitation of the method. 



result of a symmetry under an interchange of initial and final states in (^, 
and it is known as a part of Luke's theorem 0. The same symmetry also 
restricts the form of the 0{l/m?Q) terms. 

The matching factor, defined so that the identity Vo = ^vVq holds 

for matrix elements, is an ultraviolet- and infrared-finite function of mc/rrib- 
Through one-loop perturbation theory, 

r]ycb = 1 + ^Cp-— In 2 . (7) 

Att Krub — rric rric J 



The two- loop coefficient is also available p2 |. 

The vector current defined with lattice fermion fields has properties sim- 
ilar to V^ . There is a normalization factor Zy^ defined so that Vo = 
ZvqVq^^ holds for matrix elements. The factor Zvq depends strongly on the 
(lattice) quark masses arric and arrib [0, and its one-loop corrections are 



large. In the past, such uncertainties in the normalization prevented a calcu- 
lation of h^~^^{l) with the sought-after accuracy. One can, however, capture 
most of the normalization nonperturbatively by writing, with explicit flavor 
indices, 



/Zv--ZytbPycb. (8) 

In our ratio (|l]) the flavor-diagonal factors cancel, so our method avoids the 
major normalization uncertainties. 

The remaining radiative correction pyct depends on the ratio of quark 
masses and the lattice spacing. In the continuum limit, amc — >• and ami, -^ 
with mc/nib flxed, 

Pvo ^ 1, (9) 

by construction. In the static limit, anic -^ oo and anib -^ oo with a and 
rric/mb flxed, 

Pvo ^ Vv, (10) 

because the lattice theory strictly obeys heavy-quark symmetries. In numer- 
ical work one is somewhere in between, but the limits imply that py^ is never 
far from unity. Two of us have computed py^ at one loop in perturbation 
theory |^, verifying explicitly that the radiative correction is small. 
Similarly, the ratio (H) is described by the expansion 

1 - rTTT = 1 -"" + <'-'(---) - <=-' f^ - ^) + 0(™-), (11) 

n,+ (l) Krric mbJ \m^ m^ J ^ 



where j3v is a coefficient from matching the currents in (|^) to HQET. Like rjy, 
it is an ultraviolet- and infrared-finite function of rric/mh, and 

o r,r^ ^s ( '^f^biric , rrib rrib + mA ,„. 

fjycb = -i^F-r- 7 r^ In (i^j 

Air \[mi, — mc) rric nib — mc) 

at leading order. 

The ratio (p|) again captures nonperturbatively most of the renormaliza- 
tion of the lattice currents, apart from a factor pycb to compensate for the 

i 

difference between the radiative corrections with a fixed lattice cutoff and 
with no ultraviolet cutoff. In the continuum limit py- — > 1, and in the static 
limit py- —>■ 1 — j3y. Again, explicit calculation verifies that the one loop 
contribution remains small between the limits. 

In the rest of this paper, we do not write the matching factors py^ when 
there is no risk of confusion. In the final result, on the other hand, they are 
included. 



4 Lattice QCD and heavy quark symmetry 

In Ref. [|l7l, it was shown that the usual action for light quarks [l^ can be 



analyzed in terms of the operators of HQET. Therefore, it can be used as the 
basis of a systematic treatment of heavy quarks on the lattice, even when 
the quark mass in lattice units, arriQ, is not especially small. The key is to 
adjust the couplings in the lattice action so that operators are normalized 
as they are in HQET. When amq < 1, as is the case for charmed quarks 
at the smaller lattice spacings in common use, this is essentially automatic, 
because the higher order terms of the heavy quark expansion come from the 
Dirac term of the lattice action, as in continuum QCD. When amq > 1, as 
is the case for bottom quarks, one can apply the formalism of HQET to the 
lattice theory to obtain the normalization conditions, as sketched below. In 
either case, the kinetic energy is normalized nonperturbatively by tuning the 
quark mass according to some physical condition. Other operators are often 
normalized perturbatively as an initial approximation but ultimately may be 
normalized nonperturbatively. 

In the numerical calculations presented here, we use an action introduced 



by Sheikholeslami and Wohlert |16 



S = Y1 ^x V'x - Y. I^M^xyi^i + -CSW Y. ^^M^f^>^Ft^M^ (13) 

xj x,yj xj 



where the index / runs over heavy and hght flavors. The hopping parame- 
ter Kf is related to the bare quark mass, 

_ 1 1 

where K^rit is the value of k needed to make a quark massless. The flavor- 
independent matrix Mr^y vanishes except when y = x±fia, for some spacetime 
direction /i. The kinetic energy arises from this term. The gluons' field 
strength F^i, is defined on a set of paths shaped liked a four-leaf clover, so S 
is often called the "clover" action. With csw = one has the Wilson action. 

For the light quark the clover coupling csw can be chosen so that there 
are lattice artifacts of order a'^A^Qj^. In our numerical work we take an 
approximation to the optimal value, leaving an artifact of order agaAgcD- 

For heavy quarks, the clover action (|1^) has the same heavy-quark spin 
and flavor symmetries as continuum QCD, even at nonzero lattice spacing. 
Consequently, we can use the machinery of HQET to characterize the lattice 
theory. The same operators as in continuum QCD appear, but the coefficients 
can differ. Through first order in l/rriQ there are three operators in the heavy 
quark effective Hamiltonian, 

where h is a. heavy quark field, and the coefficients mi, l/m2, and l/ms 

depend on the bare mass and the gauge coupling. Because the lattice breaks 

relativistic invariance, the three "masses" are not necessarily equal, except 

as arriQ —^ 0. 

At tree level, the rest mass ami = log(l + amo), and the (inverse) kinetic 

mass 

12 1 

= 7 7 + • 16) 

am2 amQ[2 + amQ) 1 + arriQ 

The first term can be traced to the Dirac term of the lattice action, and the 
second to the Wilson term. The one-loop corrections to ami and am,2 are 
also available [^]. The chromomagnetic mass tub is considered below. 



In the heavy quark effective theory, the rest mass term mihh commutes 
with the rest of the Hamiltonian and, thus, decouples from the dynamics. 
As with decay constants [^], one can derive the expansions like @ and ([iTD 
within the lattice theory, and the rest mass disappears from physical am- 
plitudes [BH]. On the other hand, adjusting the bare quark mass so that 



^^2 = f^Q is the way to normalize the kinetic operator hT)'^h/2m2 correctly. 
This normalization can be implemented nonperturbatively by demanding 
that the energy of a hadron have the correct momentum dependence. In our 
numerical work we use the B and D mesons for this purpose. Furthermore, 
one can correctly normalize the chromomagnetic operator /iS ■ 'Qh/2inB by 
adjusting the clover coupling csw? as a function of the gauge coupling, so 
that TTiB = rn2. For example, at tree level the desired adjustment is csw = 1- 
In our numerical work, we choose csw in a way that sums up tadpole dia- 
grams, which dominate perturbation theory. This amounts to normalizing 
the chromomagnetic operator perturbatively. 

In summary, we adjust the bare mass arriQ and clover coupling csw so 
that the leading effects of the heavy-quark expansion are correctly accounted 



for IJT^. Previous work in the literature chose instead to adjust the bare mass 
until mi = rriQ, which introduces an unnecessarily large error J^ proportional 
to 1 — mi(mo)/m2(mo). 

Under renormalization the heavy quark kinetic energy can mix with the 
rest mass term in a power divergent way. Because the lattice action used here 
contains both, the rest mass fully absorbs the power divergence. A related 
problem is the ambiguity owing to renormalons |^, which appears in some 
quantities in HQET or nonrelativistic QCD (NRQCD). It is irrelevant to our 
work, because we calculate physical quantities, namely the masses of the B 
and D mesons and decay amplitude for B -^ DW. 

To complete the correspondence of the lattice theory to HQET we must 
consider the vector current. At order l/rriQ of HQET 

where the coefficient l/m^ depends on the current employed. The heavy- 
heavy current on the lattice is constructed by defining a rotated field 0,|25|, 



^f = J2^j 1 + adiiamof, gl)^ ■ D ^^, (18) 



where iIj is the quark field in the hopping-parameter form of the action (|T3D. 
Then the lattice vector current 

F/ = ^^7m^' (19) 

^To mitigate this error, these calculations are often carried out at artificially small 
quark masses. Ensuing extrapolations to larger masses contaminate lower orders in the 
(physical) l/mg expansion with higher orders. 



and V^^ = ZycbVI^^. Both Zyfg and d{ depend on the gauge couphng, the 
masses, and (at higher orders) on the Dirac matrix in (|l^. They are ad- 
justed so that the normahzation and momentum dependence of matrix ele- 
ments matches the continuum, respectively. In particular, at tree level the 
coefficient in (^) is 

1 2(l + amo) _2^^^ ^^0) 



arrij, amo(2 + amo) 

and the condition m^i = 1712 prescribes a condition on di [P^,|23|. 

From the properties of the operators under heavy-quark symmetry, it 
follows that the l/m2 and I/itlb terms in (|15|) could give a contribution 
to /i+(l), but not to /i-(l) 0. On the one hand, these contributions to /i+(l) 
must be symmetric under interchange of the initial and final states, but, on 
the other hand, they must vanish when the initial and final quark masses are 
the same. Consequently, there can be no contributions linear in either l/m2 
or l/ruB- Our definition of ^+(1) enjoys this property, by construction, 
because (|l]) manifestly respects the interchange symmetry. 

Similarly, the l/ms terms in (|T^ give a contribution only to h^{l). It 
must be anti-symmetric under interchange of the initial and final states and 
must vanish when the initial and final quark masses are the same. Our 
definition of h_{l), again by construction, ensures that only the combination 
l/^^sc ~ l/'^sfe appears. This feature is taken into account in Sec. p. 

In (^ and ( pT]) we seek contributions of order l/rriQ. These come from 
double insertions of the l/rriQ terms in ([T5|) and (|l3), and from I/ttlq terms 
implied by the ellipses. Remarkably, the latter cancel when h+{l) and h-{l) 
are defined by the double ratios (||) and (H) ||2^. This is easy to understand 
if one starts with the matrix elements. The l/ml and l/m^ corrections to 
the action arise from the initial or final state only. To this order, one can 
factorize them. They drop out of the double ratios, because the numerator 
and denominator of ([I|) or (|^) contain the same number of B and D factors. 
The same applies to 1/ml and l/m^ corrections to the current. There may be 
a nonfactorizable correction to the current with coefficient C{mc, mb)/mcmb, 
where the function C is unknown, except that in perturbation theory it starts 
at one loop and that C{m,m) = 0. 

In the long run, one would like to pick up terms of order l/rriQ and 
higher. Because the bottom quark is so heavy, these are dominated by the 
l/?7i" terms. With the normalization conditions outlined here |jT^, these 



come automatically from the Dirac term, as in continuum QCD. In future 
work, at smaller lattice spacings, the Dirac term will dominate, generating 
contributions to all orders in l/m-c. 

5 Lattice details 

Our numerical data are obtained in the quenched approximation on a 12^ x 24 
lattice with the plaquette gluon action at /3 = Q/ g^ = 5.7. We take a mean- 
field-improved [0 value of the clover coupling, which on this lattice is csw = 
1.57. Out of 300 configurations generated for our previous work [18|, we use 
200 configurations. We usually define the inverse lattice spacing through 
the charmonium IS-IP splitting, finding a~^(lS-lP) = I.I6I3 GeV. For 
comparison, with the kaon decay constant a~^{fK) = l-Oltl GeV, and the 
difference is thought to be part of the error of quenching. Because the form 
factors are dimensionless, the lattice spacing affects them only indirectly, 
through the adjustment of the quark masses. 

To investigate the heavy quark mass dependence of the form factors we 
take Kh = 0.062, 0.089, 0.100, 0.110, 0.119 and 0.125, and consider several 
combinations for the heavy quarks in the initial and final states. The mass 
of the spectator light quark is usually taken to be close to that of the strange 
quark, for which ki = 0.1405. We examine the effect of chiral extrapolation 
using four ki values, 0.1405, 0.1410, 0.1415, and 0.1419, for various combi- 
nations of the initial and final heavy quark masses Hh = 0.089, 0.110, and 
0.119. The critical hopping parameter is Kcth = 0.14327l|. 

For the computation of the matrix element^ {D{p')\V^\B{p)) we calculate 
the three point correlation function 

C^^-^(t, p', p) = 5: e-(P-P')-^e-P'^(D(0, 0)V,{t, y)B\T/2, x)) (21) 

with Vn from ([l^) and p = 0. The light quark propagator is solved with a 
source at time slice 0, and we place the interpolating field for B at T/2, where 
we use the source method. The interpolating fields B and D are constructed 
with the IS state smeared source as in Ref. |l^. The spatial momentum p' 
carried by the final state is taken to be (0,0,0), (1,0,0), (1,1,0), (1,1,1) and 

■^For simplicity we use "i?" instead of "B" to indicate the {bq) meson, and we use "B" 
or "13" for any values of the heavy quark masses. 



10 



(2,0,0) in units of 27r/L, where L is the physical size of the box; in our case, 
L = 12a. 

The numerical results presented below are obtained from uncorrelated fits 
to ratios of these three-point functions. The statistical errors are estimated 
with the jackknife method. For a subset of the data we have repeated the 
analysis with correlated fits and the bootstrap method. We find no statisti- 
cally significant difference. 

In much of the numerical work presented in this paper, we set the coef- 
ficients di of the rotation (|18D to zero. From the discussion following (|17D 
the dependence on di enters directly through l/m^, and indirectly by chang- 
ing py . On the scattering matrix elements of the spatial current Vi, this 
should make a small (^ 10 % or so) effect. On the temporal current Vq, the 
effect should be tiny. Both expectations are checked at representative choices 
of the heavy quark masses, and the uncertainty introduced into the spatial 
current is propagated to the final result. 



6 Calculation of |/i+(l)| 

The form factor |/i+(il')| at zero recoil is obtained directly from the three- 
point correlation functions (pij), setting all three momentum to be zero. We 
define a ratioQ 

W - cDVoD^i^ 0, 0)C^^oB(t, 0, 0) ' ^ ^ 

in which the exponential dependence on t associated with the ground state 
masses cancels between the numerator and denominator. When the current 
and two interpolating fields are separated far enough from each other, the 
contribution of the ground state dominates and 



i?^^^(t) 



(D(0)|K)|g(0))(i?(0)|K)|D(0)) 
{DiO)\Vo\DiO)){BiO)\Vo\BiO)) 

I + V ; + V ;i ^ |^B^D(^^^||2^ ^23) 



|/i$^^(l)/i 



B^Bl 



^Mandula and Ogilvie llfl] used a similar ratio, with nonzero velocity transfer, to study 
the w dependence of the Isgur-Wise function, which is the infinite mass limit of ft,-|_(w). 
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1.02 




t 

Figure 1: R^^^{t) as a function of t. The heavy quark hopping parameter 
for the initial and final mesons are {nb, i^c) = (0.089,0.110) (diamonds), and 
(0.089,0.119) (circles). The light quark corresponds to the strange quark, 
K; = 0.1405. The solid lines represent a constant fit with 4 < t < 8. 



suppressing radiative corrections. Here we use the definition (^) and the unit 
normalization of |/i+(l)| in the equal mass case. Thus, we expect R^^^ to be 
constant as a function of t, and its value represents the form factor squared. 

In Fig. |l| we plot the ratio R^^^{t) for two representative combinations 
of mass parameters. We observe a nice plateau extending over about five 
time slices, and our fit over the interval 4 < t < 8 is shown by the solid line. 

To see if the plateau is stable under the change of the position of the 
interpolating field, we repeat the calculation changing the time ts of the B- 
meson interpolating field. The results with ts = 10 and 8 are shown in Fig. 
together with the one with ts = T/2 = 12. We observe that the plateau is 
very stable and conclude that the extraction of the ground state is reliable. 
In the following analysis we use the result with ts = T/2, and the numerical 
data for each Kh are given in Table ^ 
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Figure 2: Check of the plateau in R^~*^{t) by varying the time shce ts of 
the B meson interpolating field. Open diamonds, crosses and solid diamonds 
correspond to the results with ts = 12, 10, and 8, respectively. The heavy 
quark hopping parameters are {k^^Kc) = (0.089,0.110), and ki = 0.1405. 



Table 1: Numerical data for i?^~*^, which corresponds to |/;,+ (l)p, at k,i = 
0.1405. Rows (columns) are labeled by the value of Kh in the initial (final) 
state. Combinations without data have not been calculated in this work. 
The diagonal elements are 1 by construction. 



I^h 


0.062 


0.089 


0.100 


0.110 


0.119 


0.125 


0.062 


1 


0.989(07) 


0.979(12) 






0.947(24) 


0.089 


0.989(07) 


1 


0.998(01) 


0.993(02) 


0.986(05) 


0.983(07) 


0.100 


0.979(12) 


0.998(01) 


1 






0.992(04) 


0.110 




0.993(02) 




1 


0.999(01) 




0.119 




0.986(05) 




0.999(01) 


1 




0.125 


0.947(24) 


0.983(07) 


0.992(04) 
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Figure 3: Chiral extrapolation of |/i+(l)p. The heavy quark hopping param- 
eters for the initial and final mesons are {nb, ^c) = (0.089, 0.110) (diamonds), 
and (0.089,0.119) (circles). 



We examine the chiral limit by computing with four values of the light 
quark mass (^^, roughly in the range 171^/2 < TUq < mg. Figure ^ shows 
that the aniq dependence of |/i+(l)p, for two combinations of (/tb, Kc), is very 
slight. A linear fit in amq gives a slope consistent with zero, and the value 
in the chiral limit is still consistent with that at the finite light quark mass. 
With our present statistics, we cannot study the dependence on the light and 
heavy quark masses simultaneously. Instead we take from Fig. ^ two lessons: 
the dependence on the light quark mass is insignificant, but the (statistical) 
uncertainty increases, by a factor of two, in the chiral limit. 

A small, but non-analytic, dependence on 771^^ is expected from chiral 
perturbation theory [^, ^. Such effects may be different in the quenched 
approximation. If so, the difference should be counted as part of the error of 
the quenched approximation. 
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7 Heavy quark mass dependence of |/i+(l)| 

In the heavy quark hmit of QCD, the heavy quark mass dependence of |/i+(l)| 
can be described with a l/mq expansion. Using a symmetry of its defini- 
tion d^) under the exchange of the initial and final states and the normaliza- 
tion in the limit of degenerate heavy quark mass, the form of the 1/mi, and 
l/rric expansion is restricted to be 



|ft+(l)| = l-cV"' +<;!;"— + — +0{l/ml 

(24) 
suppressing the radiative correction rjv- The term O^I/uiq) denotes all pos- 
sible combinations of l/rric and 1/mf, with total mass dimension —4. The ab- 
sence of terms of order l/rriQ is implied by Luke's theorem 0, but in this par- 
ticular case it can be understood as a result of the symmetry l/rric ^->- l/rrib. 
If we take the radiative corrections into account, the data presented in 
the last section correspond to |/i+(l)|/pvb- To use the right-hand side of (p4D, 
on the other hand, we must multiply by with pvo/vv to obtain \h+{l)\/r]v- 
At /5 = 5.7 and our choices of quark masses we find, at one loop, that pvo/Vv 
is very nearly 1, so that we do not need to carry out this conversion.^ 



In the lattice theory, the masses entering (|2^ are 7712, mB, and 777,3, as 
explained in Sec. |[ In particular, if one follows Refs. [^, ^ to see how 
higher-dimension tree-level operators affect the matrix elements, one sees 
that the l/'77i^ and 1/ml corrections to the action and current do not af- 
fect R^^^ |26|. 

We study the relation (|^) with several combinations of the initial and 
final heavy quark masses. We require a relation between the hopping param- 
eters, which are inputs to the numerical calculation, and the quark masses. 
To simplify the analysis, we set 7773 = 777^ = 7772 and estimate the kinetic 



quark mass by applying tadpole improvement [T^ to include the dominant 
tadpole contribution to the perturbation series. The tadpole-improved ki- 
netic mass is given by substituting arfiQ = amo/uo for a?77o on the right-hand 



side of (16), with the mean link variable Mq = l/S/tcrit- We do not bother with 



the one-loop correction to 7^72 |^3 , because it is smaller than the uncertainty 



from a. This way of parametrizing the quark masses is for interpolating only; 
when reconstituting the physical result, the hopping parameters Hc and Kb 



^This is an accident at our choice of lattice spacing. For smaller lattice spacings, this 
would not be so. See Ref. fE3| for details. 
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Figure 4: XjarHc dependence of |/i+(l)|. The initial heavy quark mass is 
fixed at K}y = 0.089, which corresponds to 1/ami, = 0.475. The hght quark 
corresponds to the strange quark, ki = 0.1405. 



are chosen nonperturbatively from the masses of the D and B mesons. 

Figure || shows the l/arric dependence of |/i+(l)|. The initial heavy quark 
mass is set to be l/arrib = 0.475 (^6 = 0.089), and we vary l/anic between 0.2 
and 2.0. (Here we misuse the meaning of subscript 6 or c to indicate the 
initial or final state heavy quark, respectively.) At l/arric = l/arrii, the form 
factor becomes exactly one by construction, and the deviation from unity 
increases as l/arric moves away from l/arrii,. The statistical error grows as 
the difference of heavy quark masses increases. When one approaches the 
static limit the signal becomes much noisier, as in many other Monte Carlo 
calculations with heavy-light mesons. In our case, the statistical error of the 
point with heaviest arric is very large. 

To see the mass dependence more clearly, we rewrite the relation (p4D as 



|/l+(l)| _ J2) J3) M , 1 



., <^-<M^ + ^ , (25) 
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Figure 5: [1 — |/i+(l)|] /A^ vs. l/arric + l/arrih. The dotted vertical line 
indicates the physical value of l/arric + l/arrib. The hght quark corresponds 
to the strange quark, ni = 0.1405. 



where A = l/arric — l/aruf,. The left-hand side is plotted in Fig. |[ The 
data exhibit a very good linear dependence on A, except in the heavy mass 
regime, where the error grows rapidly. Fitting all data linearly, we obtain 
c+ = 0.029(11) and c+ = 0.011(4). In physical units, and absorbing factors 
of a into the coefficients, these coefficients have a size typical of the QCD 
scale: c^+^ = [0.20(4) GeV]^ and cf = [0.26(3) GeV]^. 

The dotted line marking the physical value of l/arric + l/arrib shows that 
we are, in effect, using (^) as an Ansatz for interpolation. Although the 
coefficients are interesting in their own right, we caution the reader that the 
values extracted from the fit are highly correlated, and we have not made a 
full analysis of the errors on them. Below we prefer to give /i+(l), evaluated 
at physical values of the masses, as the principal result of this section. 

We have checked the influence of the rotation by repeating the calcula- 
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tions with di set to 



di 



1 



1 



(26) 



2 + arfiQ 2(1 + amo) 

at several of the heavy-quark mass combinations. This is the correctly tuned 
value at (mean- field improved) tree level Jl^. The primary effect of vary- 
ing di is through the combination (l/msc — l/msb)^ |2^. A secondary effect is 
to modify the radiative corrections factor pvo- After interpolating the masses 
to the physical point, the change on h^{l) is -1-0.00013, which is much smaller 
than several other uncertainties. Owing to this, our central value for /i+(l) 
can come safely from data with di = 0, the only value of di for which pvb is 
already available. 



8 Calculation of /i-(l) 

To obtain h-{w), it is necessary to consider nonzero recoil momentum. From 
the definition of the form factors (|^), the matrix elements of the spatial and 
temporal vector current for the nonzero recoil final state -D(p') read 



{D{p')\Vi\B{0)) = ./T?^^^ 



R^D, 



W] 



h 



B^D, 



W] 



(27) 



B^D/ 



{D{p')\Vo\B{0)) = ,/^^^^^ hX^''{w){l + w) + h^--''{w){l-w) , (28) 



where w = v ■ v' = a/T+v^, and v' = p' /rriD. 

On the lattice we start by computing the ratio of correlation functions 

%/yol*'PJ-^WoB(t,p/,0)- 
In the limit of well-separated currents, the time dependence fiattens, 

{Dip'miBm 



(29) 



r>B^D 

^Vi/Vc 



'oit^P') 



{D{p')\Vo\Bm 

f^B^D 



(30) 



2 



w] 



h 



B^Di 



W] 




h 



B^D 



h 



B^D 




w-l) 



The last step holds for small v'^ and suppresses radiative corrections. Because 
the velocity inherits statistical uncertainties from the D's kinetic mass, it is 
further useful to define a double ratio 



R 



V,/Vo 



'(t,p') = i?g7^?(t,p')/i?g74;'(t,p') 
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(31) 



Then, for large time separations, 



,(B^D)/iD^D),, ,, {D{p'm\B{0)) {Dip')\Vo\Dm 

{D{p')\Vo\Bm {D{p')\V,\D{0)) 



Vi/Vo V'Jr/ /r^/'„/\lTAlr:>/'n\\/n/„/\lT/ln/r»\\ v / 



1- 



h^-^^'ivj) 



2hx [w) 



The final expression is simplified using the property h^^^{w) = 0. Provided 
that \h^~^^{w)\ is obtained sufficiently precisely in the previous sections, the 
relation p^ ) can be used to extract h^^^(w). The part proportional to w — 1 
gives only a small contribution, since the coefficient h^{;w)/hj^{w) is itself a 
small quantity of order (m^ — mD)/{jnB + mr)). 

Figure ^ shows the t dependence of the ratio -Ry. a/q ^ (^, p') for final 
state momenta Lp'/27r = (1,0,0) (circles) and (2,0,0) (squares). Filled 
symbols represent the b ^ c transition, while open symbols correspond to 
the reverse c ^ b transition. The plateau is reached around t = 4, so that 
we can fit in the interval 4 < t < 8, as with |/i+(l)p. The fit results for the 
momentum (1,0,0) are given by the solid lines. 

Up to the small contribution of order w — 1, this ratio gives the com- 
bination 1 — h^~^^{w)/h^~^^{w), in which h^~^^{w) is almost equal to 1. 
Looking at the solid symbols in Fig. |^, h^~^^{w) is roughly —0.1 and is al- 
most independent of the final state momentum. Since h_^{w) changes its sign 
under the exchange of initial and final states, it is consistent that the open 
symbols, which correspond to the transition D —y B, appear below one. 

To obtain the value of h^~^^ (w) / h^~^^ (w) at the zero-recoil limit, we 
extrapolate the plateau values of -Ry./yg ''^ ^ ' for p' — > 0. The small piece 
of order w — 1 vanishes in this limit as well as the possible w dependence of 
form factors, so we obtain h^^^{l)/h^^^{l) without further approximation. 
Figure ^ shows the extrapolation for the same mass values as in Fig. p. There 
is no significant dependence on (ap')^. Thus, we simply apply a linear form 
to fit the data, shown in the figure. The numerical data in the zero-recoil 
limit are given in Table ^ 

The chiral extrapolation of 1 — /i_(l)//t+(l) is shown in Fig. §, for the 
combinations (k?,, kJ = (0.089,0.119) and (0.119,0.089). As in the case of 
|/i4.(l)P, the dependence on aniq is insignificant, but the (statistical) uncer- 
tainty increases, by a factor of two. 



19 



1.4 



1.2 



1.0 



DC 0.8 



; 



^ f4^U^ ^ 



4 ^Sr^ 



4r^^ 



[] 



Cl) 



0.6 



12 



Figure 6: i? 



{B^D)/{D^D) 
V^/Vo 



(t) for the final state momentum (1,0,0) (circles) 



and (2, 0, 0) (squares). The heavy quark hopping parameter for the initial and 
final mesons are {Kb, Kc) = (0.089,0.119) (solid symbols) and (0.119,0.089) 
(open symbols). The light quark corresponds to the strange quark, /t; = 
0.1405. The solid lines represent a constant fit for the momentum (1,0,0) 
with 4 < t < 8. 



9 Heavy quark mass dependence of /i_(l 

As with /i+(l), the heavy quark mass dependence of h^{l) can be described, 
in the heavy quark limit of QCD, with a l/rnq expansion. The form of the 
heavy quark expansion of h^~^^{l) is restricted by its anti-symmetry under 
the exchange of the initial and final states. 



hjl] 



---) 



a) _ ,(2) 



1 


+ 


1 ^ 


\ 


rric 


rrib^ 


K\ 



+ 0{l/ml 



(33) 



The meaning of the subscripts on the combinations of inverse masses is given 
below. The ratio /i_(l)//i+(l) obeys the same expansion up to the given 
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Figure 7: Extrapolation of Ry.T^ ~^ to the zero- recoil limit. The 
heavy quark hopping parameter for the initial and final mesons are (k^, Kc) = 
(0.089,0.119) (solid circles) and (0.119,0.089) (open circles). The light quark 
corresponds to the strange quark, ki = 0.1405. Note that the lattice spacing a 
is held fixed here. 



order, since the correction to the h+{l) starts at order I/itlq. 

To take radiative corrections into account, we should note that the (lat- 
tice) ratio i?y,/~j^p ^ corresponds to [l — Zi-Z/i+J/p^i- The right-hand side 
of (p3D, on the other hand, is justified in HQET when radiative corrections 
are ignored. Thus, we should multiply the data of Table ^hj pvj{l — Pv)- 
We shall not do this for two reasons. First, the one- loop contribution to py. 
is not yet available, although a calculation is in progress [^. Second, there 



is an indication from an analysis of renormalons that the series for Pv con- 
verges poorly [^. With these points in mind, we omit radiative corrections 
and employ (|33|) as an Ansatz for interpolation. 

The subscripts on the parentheses in ( |33D mean that the enclosed masses 
should be taken to be rris or m2, introduced in Sec. ^ The reasoning is 
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Table 2: Numerical data in the zero-recoil limit for R\r^r ~^ , which 
corresponds to 1 — /i_(l)//i+(l), at ki = 0.1405. Rows (columns) are la- 
beled by the value of Kh in the initial (final) state. Combinations without 
data have not been calculated in this work. The diagonal elements are 1 by 
construction. 



K.h 


0.062 


0.089 


0.100 


0.110 


0.119 


0.125 


0.062 


1 


1.067(12) 


1.093(14) 






1.181(21) 


0.089 


0.892(20) 


1 


1.033(04) 


1.063(08) 


1.095(11) 


1.121(15) 


0.100 


0.836(27) 


0.963(05) 


1 






1.092(11) 


0.110 




0.923(10) 




1 


1.034(04) 




0.119 




0.878(15) 




0.964(04) 


1 




0.125 


0.636(47) 


0.837(20) 


0.889(13) 






1 



as follows. The contribution to /i-(l) of first order in l/mq comes solely 
from the current 0, namely the l/m^ terms in {^7^. The second-order 
contribution comes mainly from the first-order contribution iterated with the 
l/mq corrections to the Hamiltonian |2^,|3^, namely the l/m2 and I/tub 
terms in (plsj) . We can take ttib = 1TI2 because, with the clover action, the 
difference affects the interpolation negligibly. Tracing the l/mq expansion in 
this way, and making use of the anti-symmetry under the exchange of initial 
and final states, leads to the heavy-quark expansion for the lattice data of 
the form given in (|3^) . 

In Fig. ^ we plot the l/arric dependence of /i_(l)//i+(l). The sohd circles 
are obtained by fixing the initial-state quark mass to be l/arrib = 0.475 and 
varying the final-state mass. The open circles are obtained by fixing the 
final-state mass and varying the initial-state mass. We can clearly observe 
the mass dependence, which makes it possible to extract the value of the 
form factor for physical masses. 

To extract the coefficients c_ and c_ we plot in Fig. |10 



R 



(B-.D)/(D-.D) 
Vi/Vo 

A, 



h^{l)/K{l) _a) (2)f 1 






arrir 



anihJ 2 



(34) 



where now A3 = 1/am^c — ^/afn^b. Here the solid symbols represent the 
results from the "heavier-to-lighter" transitions and the open symbols from 
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Figure 8: Chiral extrapolation of 1 — /i_(l)//i+(l). The heavy quark hopping 
parameters for the initial and final mesons are {nh, i^c) = (0.089, 0.119) (solid 
circles) and (0.119,0.089) (open circles). 



"lighter-to-heavier" transitions. The two sets of data are consistent with 
each other, except three points appearing well above the other points. These 
data involve the heaviest quark mass in our calculation, where the statistical 
noise is very large, and reliable fits become difficult. The data are well 
described by the linear form (0), and our results for its coefficients extracted 
with the "heavier-to-lighter" data are c'l^ = 0.212(31) and c'l^ = 0.054(11). 
In physical units, these coefficients are ci = 0.246(37) GeV and c_ = 
[0.27(3) GeV]2. 

The data presented in the figures and in Table ^ are obtained with the 
rotation parameter di = 0. One expects /i_ to be sensitive to di, because di 
is the coefficient of an operator of order v, and h- parametrizes a matrix 
element of order v. From the discussion in Sec. ^ however, one sees that di 
influences matrix elements through the mass rris. Thus, our method of fitting 
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Figure 9: 1/amc dependence of /i_(l)//i+(l). The initial heavy quark mass 
is fixed at Hf, = 0.089 (sohd circles), which corresponds to l/arrib = 0.475. 
The open circles are obtained by exchanging the initial and final states. The 
light quark corresponds to the strange quark, ki = 0.1405. 



compensates for the omitted rotation, provided we reconstitute the physical 
value oi h-{l) using the physical values of the quark masses throughout. A 
bonus of this method is that the radiative correction factor py^ will be easier 
to compute when di = 0. 

We have checked the infiuence of the rotation by repeating the calcula- 
tions with di = di, cf. (^). The primary effect of varying di is through l/m^ 
and, from ( ]20|) and (p3D, is proportional to the difference dl — d\. A secondary 
effect is to modify the radiative corrections of the lattice currents. 

With hopping parameters {k^, k^) = (0.089, 0.119), the difference dl — d\ 
nearly vanishes. Nevertheless, we find 



j^(B-*D)/(D- 



'^V~i) - R^^]^^^'^''^''\'^) = 0.0089 ± 0.0012, (35) 



where we use the bootstrap method to obtain a statistical uncertainty that 
takes correlations into account. This difference must stem almost entirely 
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Figure 10: — [/i_(l)//;,+ (l)]/A3 vs. l/arric + l/arrib. Filled (open) symbols 
represent the "heavier-to-lighter" ( "lighter-to- heavier" ) decay results. The 
solid and dashed lines are fitted results to the solid and open data points, 
respectively. The dotted vertical line indicates the physical value of l/amc + 
1/amfe. The light quark corresponds to the strange quark, ki = 0.1405. 



from a change in the radiative corrections, because the change in the heavy 
quark expansion is, fortuitously, negligible. Thus, it provides an estimate of 
the uncertainty from omitting the radiative corrections. 

Another check on the magnitude of the radiative corrections comes from 
comparing the heavier-to-lighter transition with the lighter-to-heavier. Be- 
cause the physical form factor h_ is anti-symmetric under interchange of the 
initial and final states, the incomplete anti-symmetry of Ry.r^ ~^ — 1, 
seen in Table |^, can come only from radiative corrections. Near the physical 
region, these discrepancies are 10-20 % of h_{l). With these considerations 
to guide an estimate, we take the uncertainty in h_{l) owing to unknown 
radiative corrections to range from +0.010 to —0.030. 
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10 Comparison with the QCD sum rules 

In the past, the form factors h+{l) and h-{l) have been studied with QCD 
sum rules or the non-relativistic quark model. Here we make a comparison 
of our results for c+ and c_ with estimates obtained with those techniques. 
From the zero-recoil sum rule, Shifman et al. obtain [^ 

,2 ,,2/1 1 \ 2 



F\ 



+ Y.Fl = l-^^^^-^{^-^] , (36) 



X 



4 \mc mi,/ 



where Fb^d corresponds to /i+(l) and the Fx represent contributions of 
higher excited states. The hadronic parameters /i^ and (Iq are estimated 
with other sum rules, and recent results are fx^ = 0.5(1) GeV^ and (Xq = 
0.36 GeV^ [H. The relation (|36D gives an upper bound for /i+(l). 



/.41)<1-^^V^ (---)', (37) 

2 Krric mi, J 

provided that the contributions F| of higher excited states are strictly posi- 
tive. This can be translated as a lower bound for the coefficient c+ : 

ct^>\{^^l-^^l) = {^.2Qtr^GeVf. (38) 

Our result c\_ = [0.20(4) GeV]^ is lower than the central value but still 
consistent within errors. 

In 0,|3O| the authors used the non-relativistic quark model to estimate 
the coefficient c+ . Their results scatter in a range (0.2-0.4 GeV^), strongly 
depending on the assumed shape of the quark-antiquark wave function and 
the value of the valence light quark mass. 



The form factor h_(l) has been studied with QCD sum rules p^,|35 
Applying their analysis to the heavy quark expansion (^) one finds 

c^'^ = ^[l + 5i-2(l + 52)r/(l)], (39) 

where A = m,B — m,b, the 6i are radiative corrections, and 77(1) represents a 
ratio of HQET form factors, at zero recoil. Neglecting radiative corrections, 
Neubert [^ finds r]{l) = 1/3 from a QCD sum rule. Taking A = 0.5 ± 
0.1 GeV and 61 = 62 = 0, this implies cL^^ = 0.08(2) GeV. With radiative 
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Table 3: Tree level estimate of the form factors at zero recoil, with statistical 
errors only. The light quark corresponds to the strange quark, Ki = 0.1405. 
The entries for the form factors do not reflect radiative corrections. 



ami) arric ^+(1) ^-(1) -^b->d(1) 

"44 n 0.992(3) -0.103(13) 1.041(8) 

3.9 1.0 0.991(3) -0.107(14) 1.042(8) 

3.4 0.9 0.990(4) -0.112(14) 1.043(8) 



corrections in the sum rule, Ligeti et al. find ri{\) = 0.6 ± 0.2 |^. Taking 
now 5i = 0.11 and §2 = 0.09 pf, this implies c^l^ = -0.05(10) GeV. Our 



result is significantly larger than both, but it is difficult to make a direct 
comparison. Our lattice calculation contains some of the radiative corrections 
automatically, and the remainder has not yet been calculated. When the 
lattice one-loop calculation is available, it should be possible to make a direct 
comparison. As we mentioned above, it is conceivable that these effects 
could change c_ significantly, without a great effect on the value we extract 
for /i_(l). 

11 Result for J^b^d{^) 

In the previous sections we have investigated the heavy quark mass depen- 
dence of /i+(l) and /i-(l) and obtained the coefficients in the l/mg expan- 
sions ( P^ ) and (P^l). To extract the value of JFb^d(I) we input the physical 
values of rric and m^, which we adjust to give the physical meson masses. 
At /? = 5.7 these parameters are arric = 1-0(1) and ami, = 3.9(5). The 
central value is fixed with the D and B meson masses with the lattice spac- 
ing a~^(lS-lP), and the error range refiects the uncertainty in the lattice 
spacing. 

The values of physical /i+(l) and h_{l) (without the matching factors) 
are given in Table ^ for three possible combinations of am^b and arric- Since 
the systematic errors in amb and in am^c are correlated, we consider the 
central and two limiting combinations only. The statistical errors on h+{l) 
and /i-(l) are estimated with the jackknife method, so that the resulting 
precision is better than that obtained by adding in quadrature the errors 
on coefficients c^ . In the physical amplitude JFb^£,(1), which is the linear 
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combination of /i+(l) and /i_(l) given in (^, the uncertainty from adjusting 
the quark masses largely cancels, and the value of J-'b-^d{^) is very stable. 

To obtain the physical result, we must now fold in the radiative correc- 
tion pvq, relating the lattice current to the continuum. Two of us recently 
have calculated this factor to one loop |2^, and at arrib = 3.9 and arric = 1.0 



they find pvo = 1 + O.OQGa^. The Lepage-Mackenzie scale q* for the coupling 
(^s{q*) dH has also been calculated, and at the same quark masses the result 
is q* = 4.4/a. At /? = 5.7, ay(4.4/a) = 0.168 and the correction to h+{l) 
is +0.016(3), taking the error of omitting higher orders to be 20 % of the 
one-loop correction. 

A similar one- loop calculation for pv^, which modifies h_{l), is not yet 
available. We allow, therefore, a systematic uncertainty for this effect. 

Our results for the form factors are 

h+{l) = +1.007 ±0.006 ±0.002 ±0.003, (40) 

h_{l) = -0.107 ± 0.028 ±0.004+°:°^°, (41) 

where the error estimates are as follows. The first error comes from statistics, 
after the chiral extrapolation; the second from adjusting the heavy quark 
masses; and the third error from unknown radiative corrections, two loops 
and higher for h^ and one loop and higher for h_. The chiral extrapolations, 
which are shown in Figs. ^ and §, double the statistical errors of Table |^, 
without changing the central values. 

Our main result is the value of the form factor entering the decay rate, 
at zero recoil. Inserting the physical values of the B and D meson masses 
and the results (g^) and (|l]) into (|), 

J^B^oil) = 1-058 ± 0.016 ± OmStomt (42) 

where errors are from statistics, heavy quark masses, and omitted radiative 
corrections. The last of these could be reduced substantially by calculating 
the radiative correction factor py. to one loop. 

Two sources of uncertainty have yet to be investigated carefully. They 
are the dependence on the lattice spacing and the effects of the quenched 
approximation. From our experience with fs |]TB|, ^ , we might suppose 
that these effects are a few percent and ~ 15 %, respectively. The ratios 
have been constructed so that all sources of error, including these, vanish 
for equal heavy quark masses. It is, therefore, our expectation that these 
percentages apply not to ^(1) but to J-'{1) — 1. That means that these two 
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sources of error should be under good control, just as we have found with 
the other sources of uncertainty. 

12 Conclusions 

In this paper we have shown that precise lattice calculations of the zero-recoil 
form factors /i+(l) and h-{l) are possible. The principal technical advance is 
to consider ratios of matrix elements, in which a large cancellation of statis- 
tical and systematic errors takes place. The numerical data are interpreted 



in a way mindful of heavy quark symmetry |17]. We find, therefore, that the 
dependence of the form factors on the heavy quark mass is well described by 
l/ruQ expansions, and we obtain the coefficients in the expansions. 

Our control over the heavy quark mass dependence allows us to deter- 
mine the individual form factors /i+(l) and /i_(l), as well as the physical 
combination TB^Di^)- The main results (^)-(|4^) account for most uncer- 
tainties, but not the dependence on the lattice spacing or the effect of the 
quenched approximation. Since our method is designed to yield the devi- 
ation of JFb^d(I) from one, we do not expect these qualitatively to spoil 
the quoted precision. With the proof of principle provided by this work, it 
should be possible, in the short term, to obtain !Fb^d{^) with control over 
all sources of uncertainty and an error bar that is small enough to be relevant 
to the determination of \Vcb\- 
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